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ABSTRACT 

We define a framework for determining constraints on the detection rate of fast tran- 
sient events from a population of underlying sources, with a view to incorporating beam 
shape, frequency effects, scattering effects, and detection efficiency into the metric. We 
then demonstrate a method for combining independent datasets into a single event rate 
constraint diagram, using a probabilistic approach to the limits on parameter space. We 
apply this new framework to present the latest results from the V-FASTR experiment, 
a commensal fast transients search using the Very Long Baseline Array (VLB A). In the 



Centre of Excellence for All-Sky Astrophysics (CAASTRO) 
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20 cm band, V-FASTR now has the ability to probe the regions of parameter space 
of importance for the observed Lorimer and Keane fast radio transient candidates, by 
combining the information from observations with differing bandwidths, and properly 
accounting for the source dispersion measure, VLBA antenna beam shape, experiment 
time sampling, and stochastic nature of events. We then apply the framework to com- 
bine the results of the V-FASTR and ATA Fly's Eye experiments, demonstrating their 
complementarity. Expectations for fast transients experiments for the SKA Phase I dish 
array are then computed, and the impact of large differential bandwidths is discussed. 

Subject headings: instrumentation: detectors — methods: data analysis — methods: 
observational — radio continuum: general — surveys 



Introduction 



Probing the time domain is emerging as one of the exciting goals of current and future radio 
instruments. Detection and characterization of fast transient events (those varying on time scales 
shorter than the correlator time scale in a typical imaging pipeline, and usually sub-second) is 
expected to probe some of the most energetic and dynamic astrophysical events, particularly for 
extragalactic sources. Pulsars are the most observed rapidly-varying radio source, where high time 
resolution voltage capture, and folding of multiple pulse profiles, are used to form an incoherent 
or coherent detection. In addition to the pulsar source population, which is well-studied, there 
are hopes for detection and characterization of other astrophysical populations, inclu ding other 
emission mechanisms from neutron stars (e.g., magnetars, RRATs, iKeane et al.l 120101 ) . sporadic 
emission from flare stars a nd planets, radio emis s ion from gamma ray bursts (G RBs), and coherent 
processes involving AGN ( Macquart et al. 201C : Cordes and McLaughlin 20031 ). Within the radio 
spectrum, the 20 cm waveband has been best studied, with two interesting ca ndidate extragalactic 
sources observed with the multibeam receiver of the Parkes radiotelescope (jLorimer et al.l 120071 : 
Keane et al.ll201ll ). These first results, and the source localization limitations of single-element 
systems, have prompted keen interest in this field. 



There have been a variety of fast transients experir nents undertaken in t he past thirty 



2009 



2011 



Stappers et al 



Burke-Spolao: 



2011 



: lBhatl 



2011 



et al.ll2011 



with a concentration of resources in the past d ecade (iDeneva et al 
Keane et allboiol 



2009; Siemion et al 



, arid plan s for future experiments (iLonsdale et al 



years 



2012 



Macquart et al.ll20ld ) . These experiments have operated in a 



range of modes, yielding constraints probing different parts of the astrophysical parameter space, as 
well as providing alternative strategies for detection. The modes include different instrument types 
(single-element and multiple-element), experimental set-up (frequency, bandwidth, channel width, 
temporal resolution), and different signal capture strategies (incoherent addition of power from 
each element, coherent addition of voltages, fly's eye pointing of elements to different fields). These 
modes aim to balance field-of-view with sensitivity, and allow optimal radio frequency interference 
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(RFI) mitigation techniques for each experiment. 



V-FASTR, The VLB A Fast Radio Transients Experiment (jWavth et al. 



2011 



2011 



Thompson et al 



Wavth et al.ll2012l ). is a commensal fast transient experiment, using the ten 25 ra dish an- 



tennas of the North American Very Long Basehne Array (VLBA) network (jNapier et alJll994l ). A 
full experiment and system description is available in these publications, and we present a brief 
description here. V-FASTR operates continuously, searching for fast transients in real time com- 
mensally with regular VLBA data acquisition. Since fields are selected based on the science goals 
of high angular resolution astronomy rather than transient science, V-FASTR dataset consists of 
observations of a range of targets with different receiver setups and system bandwidths. The re- 
sulting selection effects and biases are discussed by Wayth et al. (2012). The ten VLBA antennas 
are separated by a maximum baseline of ~8000 km, offering exceptionally high angular resolution 
when used in cross-correlation mode. Li contrast, V-FASTR incoherently adds the auto-correlated 
power from each antenna. The field of view is the same as that of an individual element, with 
sensitivity scaling as V^antS, where S is the sensitivity for an individual element and Nant is the 
number of available antennas. 

The distributed nature of the array means that interfering sources are not correlated between 
antennas, enabling excellent RFI rejection. The V-FASTR p i pelin e o ptimizes the r emova l of RFI 



using custom signal processing techniques (jThompson et al.l l201ll ) . IWayth et al.l (120121 ) present 



results for a subset of V-FASTR data (only data with 64 MHz bandwidth collected before 2012 
May). Here, we update their results for a wider range of bandwidths and more recent data. We 
also develop a framework that allows improved fast transient event rate constraints, by combining 
observations with different bandwidths. 

Apart from phenomena arising from neutron stars (pulsars and RRATs), most fast transients 
experiments have been unable to detect populations of sources. Typically, after completion of an 
experiment, an event rate constraint plot is generated, which seeks to outline the region of flux 
density versus event rate density that has been excluded by the absence of astrophysical detections. 
The transient event space is inherently multi-dimensional; in addition to source and propagation 
parameters (e.g., luminosity, distance, dispersion, time scale, sky temperature spatial fluctuations), 
there are the effects of the measurement system, which include instrument parameters (field-of- 
view, sensitivity, bandwidth, frequency), and detection parameters (incoherent versus coherent, 
single versus multiple element, sampling time scale, dispersion measure trials). The latter can 
have a significant impact on the ability to detect and confirm events, particularly with reference to 
excising radio frequency interference (RFI). 

Up to this point, there has been no clear framework for incorporating data from different 
experiments. As part of this, there is no existing framework for (1) detailed study of the impact 
of varying instrument sensitivity (frequency-dependent noise, beam shape), (2) incorporating the 
effects of experimental parameters (frequency differences, bandwidth, temporal sampling), or (3) 
incorporating detection performance degradation from the detection strategy (scatter broadening. 
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boxcar templates, RFI excision choices). In this paper we develop a framework to account for 
differences in experimental and source parameters, allowing a scaling of the results of an experiment 
to a common standardized quantity. Having scaled two experiments, we then demonstrate how to 
combine the results into a single constraint plot, using a probabilistic description of the ability 
of an experiment to include or exclude regions of parameter space. This framework is applied to 
the V-FASTR data accumulated to 2012, October 30, allowing the combination of observations 
with different bandwidths, and correctly incorporating the frequency-dependent beam shape of 
the VLB A antennas. We then demonstrate how to combine these data with the results of other, 
similar, surveys, yielding an event rate constraint cur ve in the 20 cm band with combined V-FASTR 
and Allen Telescope Array (ATA) Fly's Eye results (jSiemion et al.ll2012l ). Finally we compute the 
expected constraint di agram for the specified SKA Phase I dish system, and 461 hours of observation 
(jPewdnev et ahlboiol ). 



2. Impact of observation parameters on event rate constraints 

2.1. Bandwidth and beam shape 

We begin by defining the signal power and noise power in a time series, which has been de- 
dispersed at the correct dispersion measure. We assume the total power is the sum of the power 
in each channel, normalized by the number of channels (average power). We will initially assume 
that the source is located at the beam centre, and there is no frequency dependence to the signal 
or noise samples: 



™ 1=1 



N 



(2) 



where the sum extends over the A^ch spectral channels, and Si is the noise uncertainty in each power 
sample, in each channel. For the incoherent mean of power from multiple antennas, the signal is 
unchanged, and the noise reduced. Incorporating frequency-dependent signal and noise, extending 
to a continuous frequency coverage, and incorporating a frequency-dependent beam, B{u) yields: 

.r M A E / n^)B3{y)d^ (3) 



P. 



NchNa.nt^f^\ ~l JbW 



/ Sl,^,{v)dv, (4) 



where the integral extends over bandwidth, BW ^ the system noise, S'j^sys('^) is the noise in a 
frequency channel at frequency, i/, for antenna j, and A'ant and Az/ are the number of antennas 
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and channel spectral resolution, respectively. We have allowed for the possibility of non-identical 
antennas. Herein we assume identical antennas, with the understanding that the full expression 
can be used in the general case. 

For a candidate to be considered a detection, the signal power needs to exceed the threshold, 
given by the noise power multiplied by some signal-to-noise ratio value, C. A signal will be detected 
if it meets the following criterion: 



We represent the strength of the underlying astrophysical signal as a power-law in frequency, 
yielding P{i') = S'o(i^/t'o)°, where Sq is the flux density at reference frequency, vq. The inequality in 
equation ([5]) becomes an equality at the minimum source flux (at reference frequency vq), Smin = Sq. 
The foregoing formalism is straightforwardly generalized to include the angular dependence of the 
beam shape. The noise power is unchanged across the beam, but the signal power is attenuated by 
the beam response. The minimum detectable flux density at angle 9 from the beam centre, Smin{9) 
is given by: 



As expected, larger bandwidths and smaller thresholds yielder lower minimum detectable source 
fluxes. 

It is common to represent a homogeneous population of events by a single point in the 2D 
Cartesian plane, with the vertical axis representing signal magnitude and the horizontal representing 
rate per area per time. A survey that detects a single event provides an observed flux density and 
rate. With enough observations one could confidently estimate both the expected flux density and 
rate. On the other hand, a survey that does not observe any events provides some information to 
exclude parts of the event rate plane. For convenience, it is common to plot an exclusion boundary 
for a "null result" sky survey by the isocontour of parameter combinations that would have produced 
a single event in expectation, thereby defining a locus of Cartesian points with co-ordinates: 






(6) 



[ Rate per area per time, 5mm(^)] 



1 



(7) 




'mm 



where the area of an annulus of width has been incorporated. 
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2.1.1. Scatter broadening losses 
Due to multi-path propagation through the ISM, pulses are br oadened in time, with character - 



istic timescales that are heavily frequency-dependent — ICordes and McLaughlinll2003l ). 

Pulse broadening reduces the instantaneous signal strength, as the pulse energy is distributed over 
a longer timescale, leading to a degradation in the achievable signal-to-noise ratio, compared with 
the intrinsic pulse. For a square pulse of intrinsic width, W convolved with an exponential decay 
with characteristic timescale, r, the loss in SNR is given by: 

SNRr 



^^j, = Vl-/3 + /3exp(-l//3)^6, (8) 

where = r The scatter ing timescale in the Galaxy, r = t(DM,z/), is given empirically by 



(jCordes and McLaughlinll2003l ): 



log r = -3.72 + 0.411 log DM + 0.937(log DM)^ - 4.4 log z^gHz /US, (10) 

with substantial empirical scatter. The magnitude of the effect is dependent on the dispersion 
measure of the signal, DM, and the frequency, v. The degradation factor, e, is an approximation 
due to the scatter in the empirical relationship between DM, frequency and temporal broadening, 
and the limited availability of evidence for extragalactic sources. 

Appropriate application of this scaling relation requires a DM value to be assumed. One can 
incorporate this detection performance degradation into the event rate constraint plot. This scaling 
relation also requires knowledge of the intrinsic pulse width, W. The sampling time scale of the 
experiment, Ats, is a useful proxy for this quantity (for short pulses, the sampling time scale is the 
crucial scaling in the system). 



2.2. Intrinsic pulse width and temporal sampling 

An experiment that uses a temporal sampling time that differs from the intrinsic pulse width 
could suffer a loss of detection performance, due to the spreading of signal over time. If the sampling 
time scale is longer than the intrinsic pulse width (or, the observed pulse width after propagation 
through the plasma medium), the signal strength is degraded (averaged), and the pulse is less 
detectable. If the sampling time scale is shorter than the pulse width, detection performance 



^Obtained by integration of a square pulse convolved with an exponential tail, assuming a matched filter (optimal) 
detection template. The detection SNR for the matched filter and a general pulse shape, is given by: 



a 

where s{t) is described by a square pulse convolved with a decaying exponential scatter broadening function. 
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benefits from averaging over multiple timesteps to capture all of the pulse energy (e.g., boxcar 
averaging performed in V-FASTR, and other experiments). 

As a proxy for the signal contained in an event, we can form a composite quantity, 

^mmV^- -Sactual^/W^ Jy•S^/^ (11) 

where At^ is the experiment sampling time scale, and W is the intrinsic pulse width. This expression 
relates the true, unknown pulse width and actual source flux, to the sampling time scale and the 
measured flux, and provides a quantity that accounts for a range of intrinsic pulse widths (e.g., 
plotting the experimental minimum flux density and timescale provides an approximate measure 
of the true flux density and pulse width). Note that this expression differs from the oft-used 
energy-like quantity, SAt, by a square-root dependence on time. In this work, we are interested 
in determining the impact on signal-to- noise ratio of different effects, with a view to applying this 
framework to experiments with SNR thresholds that are used to define detections. The SNR scales 
as the square-root of the temporal sampling. 



3. A new combined metric 

We have formed scaling relations for incorporating a number of basic properties into an under- 
standing of our event rate constraints. We can now combine these to form a quantity that provides 
a measure of the true signal at reference frequency, uq (within the band), and beam angle, 9, for a 
given combination of experimental parameters {BW, C, Atg, At) and source parameters {a, DM): 



S^ctuA^o, 0; BW, C, a. At,, DM, At)Vw = ^- . (12) 

e^/N^t / ( - B{i,,e)di^ 
Jbw \^oJ 



By varying the DM, this yields a family of curves in the event rate diagram. A more general metric 
incorporates the detection performance degradation for detecting a high DM signal, and is given 
by: 

cJau [ Sl,{u)du 
e5actuai(i^o, e)^/W = T.a , (13) 

y/T^t / ( - ) B{u,e)dv 

Jew \^oJ 



As before, variation of the beam angle, 9, yields a locus of points in the event rate plane with 
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corresponding sky areas: 
[ Rate per area per t\me,eSmin,actna\(.'^o,0)VW] 



(14) 



27Te^AeAt 

min (9) 



,eS., 



mm, actual 



The event rate constraint curve plots the number of events per sky area per unit time against the 
minimum detectable flux density at frequency vq times the square-root sampling time scale (as a 
proxy quantity for the actual flux density, times the square-root of the intrinsic pulse width, times 
the performance degradation due to temporal smearing of the signal). 



4. Combining datasets 

Having demonstrated a scaling relation to take the parameters of any experiment and scale 
them to a common quantity (equation [T2|) , we can now explore how to compare and combine inde- 
pendent results into a single constraint curve. In the simplest case, where the minimum detectable 
flux density for two experiments is the same, and no transients have been found, one can simply 
add the FOVxtime quantities (survey volume) for each survey to find the total survey volume to 
be used in the abscissa of the event rate plot. However, if two surveys have different detection 
thresholds, it is unclear how the information can be combined. 

This issue highlights a broader problem in the formation and interpretation of event rate con- 
straint plots: the constraints are inherently probabilistic, and one cannot assign absolute boundaries 
beyond which a survey has excluded a region of parameter space. Instead, the noise in the dataset, 
and the presumably Poisson nature of the occurrence of fast transients, mean that the detection of 
an event is a single realization of an underlying statistical distribution. 

We begin by developing the probability distribution functions for the two plotted quantities 
in the event rate diagram: source flux density and number of events. We treat these functions 
separately. We aim to understand the chance of the null hypothesis: that there are no real events. 
Rejecting the null hypothesis at some level of significance indicates the probabitity that finite 
sensitivity and a threshold, or an unlucky timing of the experiment has caused us to not observe 
anything. 

Source flux density 

A source with true fiux density, S, is embedded within statistical noise. Hence, the measured value 
of the fiux density is a realization of a Gaussian-distributed random variable, with mean value, 
fi = S, and variance given by the dataset noise, o"^. The probability a source of true flux density S 
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will be detected above threshold, Ccr, is given by the cumulative distribution function (CDF): 

/•oo 

P{X>Ca) = / N{S,a'^)dx (15) 
1 1 jS-Ca\ , , 

where erf is the error function, and M^fi^o'^) denotes a Gaussian distribution with mean, ^u, and 
variance, a'^ . The cumulative probability in equation 1161 describes the exclusion regions. The cu- 
mulative probability that an event is not detected (due to noise) is the complementary function, 
1 - P{X > Ca). 
Expected event rate 

Real astrophysical events are assumed to be randomly distributed on the sky, with a given (un- 
known) mean frequency of occurrence, but random actual timing. This distribution is governed 
by the law of rare events, and follows the Poisson distribution (observation of zero fast transient 
events during an experiment does not imply that none would have been observed in an identical 
experiment). The probability of observing k events, given a mean (expected) number. A, is given 
by: 

P(.;A) = ^^^1^ (17) 

= exp-A \k=o, (18) 

implying that for k = (zero observed events), there is a non-zero probability of seeing one or more 
events in an identical experiment. In the event rate diagram, the Poisson probability distribution 
is plotted at, 

X = \/{FOV X At), (19) 

where X is the expected number density of events (compared with the observed number, typically 
plotted in an event rate diagram). If we have observed zero events, the chance that we should have 
observed at least one is given by the complementary function, 1 — Pr{k = 0; A). 



4.1. Allowed region of parameter space 

There are two possibilities when interpreting a null result experiment: (1) that no events were 
observed because there was nothing there (the null hypothesis), or (2) that fast transients are 
present, but noise or stochastic timing prevented observation. These two hypotheses have a total 
probability of unity. The probability of the null hypothesis, P{Ho), corresponds to the chance that 
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at least one source would have been detected, and is given bjo: 

P{Hq) = at least one signal is above the threshold when n events do occur (20) 
= Y^il -P{X< Cay-) x{P{k = n;X)), (21) 

n>l 

and the probability that events do occur, but were not detectable, is given by: 

1 - P{Ho) = 1 - ^(1 - PiX < CaT){P{k = n- A)). (22) 

n>l 

Figure [T] demonstrates the new event rate plot. The quantity plotted is the probability that an event 
in this part of parameter space would not be detected by our experiment (hence, low probabilities 
correspond to traditional areas of exclusion — frequent, bright events). At the beam centre, where 
the signal is not attenuated, the experiment is able to exclude regions at lower flux densities, but 
at the expense of small FOV (larger expected event rate required). 

The combined probability for multiple experiments considers the joint probability that events 
occur but are not observed in either experiment. The total probability for N independent experi- 
ments is given by: 

N 

PTot = '[l{l-Pi{Ho)). (23) 

i=l 

This expression treats the experiments independentl}{§, and creates more stringent event rate con- 
straints as more information is added. For a weak experiment (high minimum flux density, for 
example), the null hypothesis has zero probability for a large amount of parameter space, yielding 
little additional constraint when combined with a stronger experiment. Note that this expression 
may not give the same answer as taking two non-identical experiments and summing the total 
survey volume for a stronger constraint. This is because the probability description includes the 
sensitivity of an experiment, which, in general, is not the same for two experiments. 

We now have all the tools to be able to combine constraints from different experiments. It 
remains to discuss the types of experiments that are meaningfully combined, and those that are 
not. In general, any experiment could be combined with any other, and all of the differences 
incorporated into an overall framework. However, in practise our new framework only considers 
noise level, beam shape and signal differences. Results from the same instrument and different 



^Mathematically, this expression sums over possible non-zero true events, and multiplies the probability of n 
occurring by the probability that at least one would be above the threshold. 

^Independence is appropriate for the case of combining information from different pointings, or different instru- 
ments. For constructing an event rate plot for the same pointing, where different areas of the beam have different 
minimum detectable flux densities, only the new part of the field-of-view is independent (the incremental information 
added by a small area is independent of the information available from the existing area use in the plot). To con- 
struct an event rate constraint diagram across the beam, we treat each area annulus independently, and incrementally 
combine information. 
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Probability V= 1 , C = 5, cr=l 




-1.5 -1.0 -0.5 0.0 0.5 

log Expected event density (deg^ . h)~' 

Fig. 1. — Example probability event rate plot, displaying the probability that events are present, 
but not detected (due to noise or the stochastic nature of event timing). Threshold is set to C=5 
with cr = 1, At^ = 1 s, FOV.Ai=l deg.^s. 
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bandwidths can be meaningfully combined, but this framework may be too simplistic to combine 
single element with multiple element experiments. This is because it considers the data to be clean 
of RFI, justifying the use of Gaussian statistics to describe the data. In practice RFI can hurt an 
investigation in one of two ways: first, by causing false positive detections, which does not affect 
the sensitivity limit (though it may raise the burden of candidates to examine); and second, by 
artificially raising the background signal used to compute Ca, and as a result reducing the actual 
sensitivity. In general, RFI excision with single element experiments is difficult, due to a lack of 
independent information, and datasets typically contain some level of RFI. Indeed, even V-FASTR 
is subject to RFI, and occasionally RFI spikes add incoherently to form a detection. The different 
strategies employed by these experiments lead to differing RFI statistics in the data. Because our 
framework does not consider these aspects, it cannot provide a level playing field to combine all 
experiments. 

This highlights a more general discussion about the performance of a detector. In this work, 
we have focussed our attention on forming a more realistic event rate constraint diagram, with 
no reference to the design and performance of a given detection strategy (we have simply taken a 
threshold signal-to-noise ratio). In practise, the threshold, C, is chosen to balance false positive 
detections with false negatives (true astrophysical transients are not detected). The overall perfor- 
mance of a detec ti on str ategy can be quantified by exploring this balance as the threshold is varied. 



Thompson et al.l (j201ll ) explore the use of receiver operator characteristic curves (plotting false 



positive fraction versus false negative fraction) to quantify the performance of different detection 
strategies. In the V-FASTR context, antennas are combined in different ways to discriminate true 
sign als from RFI. Th e data are inherently non-Gaussian, due to the underlying RFI environment, 
and [Thompson et al.l (120111 ) incorporate this information into the data statistics to form useful 
detectors. In a similar way, the assumption of Gaussian-distributed noise used here (equation [T6]l 
could be extended to use any general distribution function (including the statistics of remaining 
RFI), and the event rate constraint curve altered to reflect thij^. With a good understanding of the 
RFI environment of each experiment, a broader class of experimental results could be combined 
meaningfully. This generalization is beyond the scope of this work. 

We use the description in equation [12] to standardize each experiment, taking into account 
performance losses due to scattering, bandwidth and sampling time scale. We then form the proba- 
bilistic description of the event rate constraint, according to equation 1221 and combine according to 
equation 1231 Note that this same formalism is used to combine results from different beam angles, 
for the same experiment. The validity of this assumption rests with the assumption that all look 
directions are equivalent. 



^V-FASTR does have subtle non-Gaussian effects due to excising large-magnitude datapoints during RFI- 
contaminated observations. However, there is not a well-defined analytical way to model the resulting sensitivity 
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5. Application: V-FASTR constraints 



The V-FASTR experiment has produced event rate constraints for ten observing wavebands, 
ranging from central frequencies of 325 MHz to 86.2 GHz. Over this frequency range, the beam 
size, SEED, and time-on-sky, vary considerably. We will first explore how the previously-used 
event rate plot is modified by the inclusion of experimental parameters, and then produce the new 
plots for four interesting wavebands. As with previous estimates, we assume that the antennas are 
identical. For simplicity, we will also assume that there is no variation of system temperature over 
the bandwidth of an observation (although we have developed the framework to account for this). 

Table [1] shows the total time-on-sky results, by receiver and bandwidth, for V-FASTR, as of 
2012 October 30. The total time is divided into three bandwidths, BW = [64 MHz, 128 MHz, 
256 MHz]. Some frequencies have substantial time at the larger bandwidth. We can use the 
framework developed here to combine the information from different frequencies. 

The VLB A antennas are 25 m dishes, with beam patterns that are well-described by Airy disks 
dohawanlbooa ): 



B{u,e) oc 



J^(u) 



cJi{'kDu sva.9 / c) 
TiDv sin 9 



(24) 



where Ji{u) is the Bessel function of the first kind. The Airy function passes through multiple nulls 
off-axis, and a source located close to a null will have a corresponding large minimum detectable 
flux densit}{§. We truncate the VLBA beam beyond the t hird Airy null (tt ~ 13.3) , and set the 



gain t o a constant value to the horizon, following the work of lDeneva et al.l (|2009l ) and lWayth et al 
Eli). 



^The beam pattern of a VLBA antenna is well approximated by the Airy function, which possesses a series of nulls 
which may be determined approximately from the asymptotic form of the Bessel function, Jm{z) for z S> \m? — 1/4|: 

Jm(2) « J— cos (2 - ^ - . (25) 

V TTZ V 4 4/ 

Thus the Airy function possesses a null at the points x = jitt, where n = 1, 2, 3, . . .. (There is no null at n = since 
limi,_>o(2Ji(a;)/a;)^ — 1.) It is evident that the spectrum of a source that is observed at an angular position 6q well 
off axis (i.e. > 1/ka) and over a large bandwidth will be subject to a large number of oscillations, due to the 
oscillations in the beam response over the wide range of k (wavenumber). We can quantify this by examining the 
number of nulls imprinted in the observed source spectrum due to oscillations in B{9). Since a null occurs every time 
kaOo changes by vr, the number of nulls for an observation between frequencies i^i and V2 is 

^ _ k2a9o - fciagp 

TV 

Note that A'^ depends only on the total bandwidth, and not on the relative bandwidth Af/u. Thus the number of 
oscillations present in the detection spectrum indicates the offset of the source from the pointing centre. However, for 
a fixed ^0 obviously the beam width is smaller at high frequency, and so for a source well off ax;is the overall amplitude 
of B{9) will be much smaller at high frequencies, rendering the source harder to detect, for a given (constant) flux 
density. 



26oa, 



[V2 



^1) 



(26) 
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We first present the latest set of results for the 20 cm receiver using the existing event rate 
fr a.mework, to demons trate the improvement in time-on-sky compared with the results presented 
in IWavth et alj (|2012l ). Figure [2] displays the constraint curves for BW = 64 MHz and different 
dispersion measures (a), and for all bandwidths and zero dispersion measure (b). We also plot 
the existing constraint for BW = 64 MHz. The curves track as expected — longer on sky yields 
stronger constraints on the event rate, while the larger bandwidths yield lower minimum detectable 
flux densities. The contributions from the three inner sidelobes are evident as cusps in the constraint 
curves. 

Figure [3] displays the probability contours for the 20 cm results from V-FASTR, plotting the 
probability a candidate in that region of parameter space would not be observed by the experiment. 
The abscissa is the expected number of events per square degree per hour, and the ordinate is the 
combined detection quantity: eS^ctnaiVW (with dimensions, Jy-^/s), formed according to equations 
[13] and [23l We also display the event rates for tw o experiments with t he Parkes radiotelescope , 



which both found a candidate fast transient event (iLorimer et al.l 120071 : iKeane et al.ll201ll . l20ld ). 

The measured DMs for these candidate events have been incorporated into the efficiency factor, e, 
and used in the plot. The structure of the isocontours is real, and corresponds to the inner sidelobes 
of the VLBA antennas. Table [2] displays the values used to plot these events. Interestingly, the 
event rate implied by the Lorimer event is on the cusp of expected detectability for the 20 cm 
receiver of V-FASTR, and future data releases will begin to challenge detections in that region of 
parameter space (assuming no V-FASTR detections). 

In figure S] we display the constraints for other interesting V-FASTR receiver bands. All of 
these wavebands have substantial time-on-sky at the higher bandwidth (256 MHz). The longer 
wavelengths sample a larger field-of-view, corresponding to a greater survey volume, and lower 
constraints for the same time-on-sky. They also have lower system temperatures (SEFD), implying 
better sensitivity constraints. In particular, the 4 cm receiver has the advantageous combination 
of long time-on-sky, low system temperature, and large field-of-view, yielding strong constraints at 
~8 GHz. 

We now use the new framework to explore the combined constraints from different surveys. 
As discussed previously, it is not reasonable to form a sensible comparison of single-element and 
multiple-element experiments, due to the different RFI mitigation strategies one must employ. 
Instead, we choose to combine the V-FASTR results with those from the ATA Fly's Eye experiment 



(jSiemion et al.ll2012l ). In 'Fly's Eye' mode, the ATA employs 30 of its 6-m dishes, each pointed in 
a different direction, leading to large instantaneous fields-of-view, but limited sensitivitjo. This is 
in constrast to V-FASTR, where the incoherent combination of dishes lead to good sensitivity, but 
limited field. Figure [5] displays the event rate constraint diagrams for the ATA alone (a), and for 
the combined ATA -|- V-FASTR experiments. The different strengths of each experiment combine 



Assuming an Airy disk for the ATA beam shape, which may not be a good approximation. Different beam shapes 
will alter the details of the cusps in the event rate diagram, but not the overall shape (which is set by the dish size). 
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Receiver 


AtTot (h) 


At64 


Ati28 


At256 




SEFD (Jy)^ 


90cm 


8.33 


8.27 


0.06 


0.00 


9.52 


2227 


50cm 


0.06 


0.00 


0.06 


0.00 


10.00 


2216 


20cm 


704.57 


539.87 


43.53 


121.17 


9.44 


296 


13cm 


197.38 


188.52 


4.27 


4.60 


9.03 


322 


6cm 


48.11 


27.99 


9.20 


10.92 


9.65 


210 


4cm 


732.99 


423.73 


215.68 


93.57 


10.13^ 


307 


2cm 


385.47 


287.36 


6.09 


92.02 


9.57 


550 


1cm 


480.57 


411.98 


43.27 


25.31 


9.83 


502 


7mm 


328.01 


253.86 


2.41 


71.74 


9.45 


1436 


3mm 


29.91 


17.92 


0.00 


11.99 


7.72 


4000 



Table 1: V-FASTR time-on-sky results for each receiver, as of October 30, 2012. The total 
time is split into time with bandwidths, BW = [64 MHz, 128 MHz, 256 MHz]. Also shown 
is the average number of antennas available with that receiver, A'^ant, and the SEFD used 
in the noise calculations. 



"Obtained from http://www.vlba.nrao.edu/astro/obstatus/current/node7.html 



''At 4 cm, additional antennas are sometimes available, when the VLB A is used for 
geodetic purposes. 



Candidate 


DM Ats ^peak AtTot.FOV 


Unit 


pc.cm"^ s Jy deg.^ h 


Lorimer et al. (2007") 
Keane et al. f2011) 


375 0.005 30 1303 
745 0.007 0.41 4264 



Table 2: Event rates from the two candidate events from the lLorimer et al.l (120071 ) and lKeane et al 
(|201ll ) experiments with the Parkes radiotelescope (note that the thirteen multibeam receivers have 
been incorporated into the survey volume metric, Aixot-FOV). We have assumed an Airy disk, and 
extended to the third null, as with the VLBA antennas. 
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(a) Variation with source DM (BW = 64 MHz). 

A = 20cm, All BW 



64MHz, 540h 
128MHz, 44h 
756MHz, 1?1h 
64MHz, 245h 




\-8 



0"^ 10"^ 10" 

Event rate (/degyhour) 



,-4 



>-2 



1 0^ 



(b) Variation with bandwidth and observing time (DM = 0), including the previous result from 
Wavth et all |2012h (black, dash-dot-dot-dot). 



Fig. 2. — Constraints from recent V-FASTR accumulated data, using the previously- used event 
rate plot, and the 20 cm receiver. The Airy disk beam has been truncated after the third sidelobe, 
with constant gain to the horizon beyond that angle, (a) Limits for observations with 64 MHz 
bandwidth (A t=540 h) and differ ent source DM. (b) Limits from each bandwidth (DM = 0), and 
the hmit from lWavth et al.l (|2012l ) (black, dash-dot-dot-dot). 
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Probability non — detection, A = 20cm 




-6 -4 -2 

log Expected event density (deg^ . h) 



Fig. 3. — Event rate constraint curve for the 20 cm receiver of V-FASTR, using the probabihty 
framework an d the new detection metric (eSar.tm ] VW) . Also plotted are the event rates for the 



Lorimer et alJ (120071 ) (asterisk) and iKeane et alJ (|201ll ) (plus) experiments, incorporating loss in 



detection performance due to dispersive temporal smearing. 
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(a) 4 cm receiver. 
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(c) 1 cm receiver. 



Probobiliiy non-detection, A = 2cm 
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(b) 2 cm receiver. 
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(d) 7 mm receiver. 



Fig. 4. — Event rate constraint curves for four interesting wavebands, with the V-FASTR experi- 
ment. 



to form a strong con strai nt across much of p arameter space. Again, we also plot the results from 
Lorimer et al.l ([20071) and iKeane et alj (|201ll ). The limited sensitivity of the ATA dataset cannot 
add significant information at these flux density levels. 



5.1. Fast transients constraints with SKA Phase I 



The system specifications for the SKA Phase I dish array outlin e an instrumeiit with 250 
dishes distributed over an area with a maximum baseline of 100 km (jPewdnev et al.ll20ld ). In 
the high-band (1-2 GHz), the description specifies a maximum bandwidth of 1000 MHz, yielding 
an instrument with su bstantially di f ferent beam shape characteristics at each end of the band. 
Following the work in IWavth et al.l (120121) . w e generate an event rate constraint curve for the 
system specifications in iDewdnev et al.l (|20ld ) in the high-band, using an incoherent addition of 
antennas, and 461 hours of observation (figure [6]). The sidelobe pattern is again ob served in the 
structure of the isocontours. The plot deviates from that shown in lWavth et al.l (|2012l ). due mostly 
to the significant fractional bandwidth, and the correct accounting for sensitivity variation with 
frequency and angle from the boresight. In addition, the SKA Phase I dish dump rate is 0.1 s (due 
to limitations imposed by the long baselines), yielding degraded performance for fast transients on 
shorter timescales. This is incorporated into the metric plotted here. Compared with the ATA, 
the SKA has substantially better sensitivity, but poorer field-of-view in incoherent addition mode. 
SKA d ishes may also be dep loyed in other fast transient configurations, and these will have different 
limits (jColeeate et ahlboid ). 



6. Conclusions 



We present the latest set of results from the V-FASTR experiment, a commensal fast transient 
experiment using the VLB A network across ten wavebands. These results are presented through 
the lens of a new framework, developed to incorporate source and experimental parameters into 
a metric that (1) captures additional information that is meaningful for characterizing the true 
sensitivity of an experiment, and (2) allows independent datasets to be combined in a rigorous 
manner. We use this framework to combine results from different bandwidth observations within 
V-FASTR, demonstrating that the 20 cm receiver o bservations are pushing into the para meter 

). We 



space defined by the Lorimer burst and Keane event (jLorimer et al 



2007 



Keane et al. 



2011 



then combine V-FASTR results at 20 cm with the ATA Fly's Eye dataset, demonstrating the 
complementary nature of the two experiments. Finally we present expectations for the SKA Phase 
I dish array instrument, producing plots that deviate from those presented previously, due to the 
correct accounting for beam shapes as a function of frequency and angle from the boresight. 
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(a) ATA Fly'e Eye constraint curve dSiemion at aLlboid V (b) Combined 20 cm limit from V-FASTR and ATA. 



Fig. 5. — Event rat e constraint curves for the ATA Fly's Eye survey, using the new framework 
(jsiemion et al.|[2012l ) (a), and combined with the V-FASTR results (b). 




Fig. 6. — Event rate constraint curves for the SKA Phase I dish array, assuming 1 GHz bandwidth 
between 1 and 2 GHz, 250 dishes in incohe rent mode, and 461 hours of observation (following 
system specifications of lDewdnev et alJl201Cll ). 



